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1.  INTRODUCTION 

Let  X be  a compact  subset  of  the  real  line  and  denote  by  C(X) 

the  class  of  all  real  valued  continuous  functions  defined  on  X.  Norm 

C(X)  with  the  uniform  norm,  l.e.  for  all  f e C(X), 

||f||  ■ max{|f(x)|:  x e X) . Let  n be  a positive  Integer  and  set 

R^(X)  m {1/p:  pell,  p(x)  > 0,  Vx  e X}  where  n denotes  the  set  of  all 
n n n 

real  algebraic  polynomials  of  degree  s n.  Note  that  R^(X)  consists  of  the 

positive  elements  of  the  usual  setAdenoted  by  Rn(X). 

In  this  paper,  we  shall  study  the  problem  of  approximation  of  positive 

functions  in  C(X)  by  elements  of  R^(X).  The  emphasis  of  this  study  is  as  follows. 

First,  we  wish  to  contrast  this  setting  with  that  of  approximation  by  elements 

of  R^(X),  m _>  1,  n _>  1,  in  C(X).  Basically,  there  is  one  major  difference, 

namely  that  existence  holds  for  this  case;  whereas,  this  is  not  true  for 

Rm(X),  m > 1,  n > 1 and  X not  an  interval.  In  addition,  when  X is  not  an 
n — — 

interval  the  proof  of  existence  is  very  long  and  tedious.  Next,  we  shall 
observe  that  the  usual  characterization  (alternation)  and  uniqueness  results 
hold  for  this  problem  using  the  standard  arguments.  Finally,  we  shall  discuss 
the  computation  of  best  approximations  from  R^(X). 

These  results  will  be  used  in  a forthcoming  paper  on  uniform  approximation 
on  [0,  ®)  with  reciprocals  of  polynomials.  See  ([2],  [3],  [4],  [8],  [14]  and  [18]) 
for  various  recent  studies  concerning  certain  aspects  of  this  problem.  Since 
we  are  interested  in  applying  these  results  to  this  setting,  we  shall  approximate 


1/f  in  what  follows. 
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2.  EXISTENCE 


As  noted  in  the  introduction,  it  is  well  known  that  best  rational 
approximants  from  R™(X) , n 2 1,  n 2 1 and  X not  an  interval  need  not 
exist.  Indeed,  let  X be  a finite  subset  of  [0,1],  such  that  0 e X 
and  card(X)  > n + m + 2.  Define  f e C(X)  ~ R™(X)  by  f(0)  = 1 and 


2kx+l 


f(x)  * 2 for  x e X ~ {0}.  Let  r^(x)  be  defined  by  r^Cx)  ” 


Then  lim  f - r. 


0 showing  that  no  best  approximation  exists  for 


f from  R (X) . 
n 


In  what  follows,  let  n be  a nonnegative  integer  and  let  X be  a 
compact  subset  of  [0,B],  B > 0 with  card(X)  2 n + 2.  Set 
K ■ {p  e nQ:  p(x)  > 0,  Vx  e X),  let  f e C(X) , f i K with  f(x)  > 0 

for  all  x e X and  set  A » inf{|  |y  - | : p e K} . We  have  the  follow- 

ing lemma. 


LEMMA  1.  Let  X,  f,  A be  as  defined  above.  Then  A > 0. 

2.  c 

Proof.  Suppose  A * 0.  Then,  setting  — ■ — , c a constant, 
q e K,  | | q | | ■ 1,  there  is  a sequence  of  polynomials  {qn>  c K such 

1 c 1 

that  12  1 1— -|  | -►  0 . Since  | c | £ 1 + — , m * min{  | f (x)  | : x e X} , 

qn  n m 

we  may  assume  that  cq  -*•  c*  2 0 and  qQ  -*■  q*  uniformly,  with  q*  e 11^ 
and  | |q*| | ■ 1.  We  claim  that  c*  > 0.  Suppose  c*  ■ 0,  and  let 


z e X be  such  that  q*(x)  > 0 (such  an  x exists  since  q*  t 0 and 
card(X)  2 n+2) . Then  we  have 
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0 lim  1 


n <-  1 . 1 0 _ lim  1 Cn  lim  i 1 1 Cn  i i ... 

f(x)  f(x)  ” q*(x)  ir*”  f(x)  qQ(x)  n-*®  f qn  " C 


is  a contradiction.  Hence  c*  > 0.  Since  | | — | | s 1 + | |-j|  | we  have 


cn  c* 

that  q*(x)  > 0 for  all  x e X.  Thus ► — -r  e K uniformly  and 

qn  q 

| |4  - -^r|  | * |-^ — | | ” 0,  hence  f K,  which  is  a contra- 

r * r q 


diction.  Therefore  A > 0. 


THEOREM  1.  Let  n b£  a nonnegative  integer  and  let  X,  K,  f and 
A be  defined  as  above . Then  there  exists  a p*  e K such  that 


llj-^ll  -a. 

£ P 

Proof.  Since  X is  compact,  f > 0 on  X,  there  exist  positive 

constants  m,  M such  that  m £ — <,  M for  all  x e X.  For  p e K 

f (x) 

we  write  * q where  q(x)  > 0 for  all  x e X,  c > 0 and  | |q| | ■ 1. 

Then,  as  in  Lemma  1,  there  exist  sequences  {cq} , {q^}  such  that  cq  > 0 
for  all  n,  q^Cx)  > 0 for  *11  * e X with  I 1^1  I * 1 satisfying 
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i)  1 + A * | |i  - f\  | 
Hn 

ii)  11m  | l)  - ^1  I -A 
n-*»  4n 

ill)  | c | S 1 + A + M. 
n 


jxM& 


*u 


By  extracting  subsequences  {c  } 

: — 


respectively,  there  exists  c*  £ 0 and 


x e X and 


1 such  that  c 


and  (q  } of  {c  } and  {q  } 
nv  n “ 

q*  e II  with  q*(x)  £ 0 for  all 

-*■  c*  and  q (x)  ■*  q*(x) 

ft. 


uniformly  in  X. 


We  now  claim  that  c*  > 0.  Indeed,  suppose  c*  - 0 and  let 
Z m (x^ x^}  <=  X,  k < n,  be  all  the  zeros  of  q*  contained  in  X. 

Now,  if  Z * 0,  then  q*(x)  > 0 for  all  x e X holds  and  we  have  for 
each  x e X 


f(x)  f(x)  q*(x) 


0 -lim  1 


s lim||i- -^||  - A. 

\ 


f(x)  Si  (X)  ^ 

V 


(1) 


~ 2 

Thus,  taking  p " — , the  inequality 


1 fYxy  - — I “ 1 - yl  * j holds  for  each  x e X.  Since  p e K, 

this  contradicts  our  assumption  that  inf  j|  |-j  - — | | j p e K>  * A > 0. 

pell  ''  p J 


Therefore,  we  assume  that  Z 1*  0.  Partition  Z into  two  subsets 


Z ■ I u J where  I ■ (x  e Z:x  Is  an  Isolated  point  of  X) , and 
J - Z ~ I.  Now,  for  x e X ~ Z we  have  that  q*(x)  > 0 so  s A 

holds  by  (1).  Also,  since  f e C(X)  and  y e J implies  that  y is  a 
limit  point  of  X ~ Z we  have,  by  continuity,  that  5 A holds. 

Now,  if  1=0,  then  our  argument  for  the  case  Z = 0 yields  our 
desired  contradiction.  F nee,  let  us  assume  that 

I • {y^  < y2  < < y^},  y £ 1.  Now,  let  a < b be  such  that  X c (a,b) 

Construct  open  intervals  (a^,8^)  as  follows  (observing  that  X~  I is 
a compact  subset  of  X):  set  a.  • max{a,  max{x:  x e X ~ I and  x < y.}} 


and  6^  ■ min{b,  min{x:  x e X~  I and  x > y^}}.  If  8^  > y^  stop  this 
process  with  the  interval  If  8^  < yy,  then  there  exists  an 

integer  i^,  1 s ij  s y-1  such  that  y^  < < +1*  In  case, 

set  a2  ■ max{x:  x e X~  I and  x < y^  +^}  and 

$2  ■ min{b,  min{x:  x e X ~ I and  x > y^  +^H.  Note  that  we  must 
have  ®2  ^ 88  8,  £ X~  I and  8-,  < y.,  Al . Once  again,  if  6,  > yti , 


1 7i,+l* 


2 V 


stop  this  process.  If  B2  < yy»  we  continue  and,  since  y is  finite, 
this  construction  must  end  after,  say  v £ y steps,  giving  v pairwise 
disjoint  open  intervals  (a^,8^> , . . . , (ay,8v)  where 

* 5 al  < ®1  S a2  < 02  5 *“  5 av  < 5 b>  n 1 1*  0 and 

n (X  ~ I)  =0  for  i ■ l,...,v.  For  convenience,  let  us  assume 

that  (ar.8r)  n I - {y^  +1 y±  > for  r - with  iQ  - 0. 





Now  we  will  construct  a new  set  of  points  I’  where  I’  n X ■ 0 as 

follows:  if  (o  ,B  ) is  such  that  a < a , $ < b and 

r*  r r*  r 

(0^,6^)  n I - {y^  +^,...,y^  ^ consists  of  an  odd  number  of  points, 

then  set  ■ -^(8r  + yi  ^ an<*  require  that  y^  el'.  Therefore,  I' 

r r r 

can  consist  of  at  most  v points.  For  convenience,  set 

1'  * {y| y^},  y ^ v,  if  I'  ? 0.  Also,  in  this  case,  note  that 

q*(x)  must  vanish  at  y and  q*(ar)  > 0,  q*(8r)  > 0 both 

r-l  1 r r r 

hold  implying  that  q*  must  have  either  one  y^,  i^  ^+1  ^ i ^ if  as  a 

zero  of  even  order  (at  least  two)  or  have  at  least  one  more  simple  zero 
in  (ar,8r>.  Thflt  is,  q*  must  have  at  least  if  - ip_^  + 1 zerosin 

(ctr,Br).  Set 


v(x) 


n (y'  - x) 
i-1  1 


I ft  0 


(2) 


N 

Next,  we  shall  construct  a sequence  of  polynomials  ^®^jj=1 


corresponding  to  each  interval  (a^.B^),  i = l,...,v.  First,  set 


/ j-i  y 

I*),  - f(yJ/N  n (y  - y ) n (y  - y.)|v(y  )| 
j ->  i=l  J 1 i-j+1  1 j 3 


(3) 


0 y 

for  J ■ l,...,y  where  n = II  = 1 and  note  that  w.  0 as  N 

1 y+1 . 3 
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Let  us  first  consider  the  interval  (a^,B^)  where 

M 

(aliSl)  n I ■ {y1,...,y1  },  1 ^ i^  s y . The  precise  form  of  0^  will 

depend  upon  the  structure  of  although,  in  all  cases,  the  poly- 

nomials  0^  will  have  certain  essential  properties.  Therefore,  we  must 
consider  cases. 

Case  1.  = a,  i^  = 2m  + 1,  m i 0.  Note  that  in  this  case  we 

must  have  B^  < b since  card(X)  2 n + 2.  Furthermore,  this  interval 
gives  no  contribution  to  I ' . Set 


0j(x) 


m 


m 


- n [(y.,,.,  - (0  ) -x]  n [(y_.  + u>,,)  - x]  where 


j=0 


2j+l  2j+l 


j-1 


2J  2j 


is 


defined  by  (3).  Now,  since  u>j  -*■  0 as  N -*•  » for  all  j and 
a » < y^  < • • • c y^  < B^,  we  can  select  an  N1  such  that  N 2 N^ 


implies  that 

4 yl  “ “l  < yi  < y2  < y2  + “2  < ***  < 

n 

Thus,  for  N 2 N^,  we  have  that  0^(x) 


yiri * yiri + "1 
is  positive  at  y 


1 

1 


-1 


0^(x)  > 0 for  x 2 B^  (which  implies  that  0^(x)  > 0 for  all  x e X) . 
Let  p(s)  denote  the  number  of  zeros  of  q*(x)  in  the  interval 
(c*a,Bs)»  s ■ 1, . . . ,$  where  we  count  a zero  of  order  p as  p zeros. 
Then  we  have  30^  s p(l)  in  (a^,6^).  Also,  setting  * y^  - for 

i odd  and  * y^  + for  i even  and  noting  that  e^  -»•  as  N -*■  ®, 

we  have,  for  Is  is  i^. 


that 


10 


, JJ  4X1  1 1 

1/NMyJ  = 1/N(-1) 1 V n (e.  - y.)  n (e.  - 

J-l  3 J-i+1  j 


?i> 


1-1 


i!i(yi  ~ i'j)1,nltl(yj ' yi)|vM 
i-1  il 

f(y1)  it  (yA  - e ) n (e  - y ) 
J-l  3 j-i+l  J 1 


Now,  using  (2)  and  (3) , we  have 

l/N0^(y  ) -+■  f(y  ) 1 n (y  - y ) n (y'  - y ) 
j-ij+1  3 j-l  3 1 


as  N 1 £ 1 i 1 


Also,  for  x e X satisfying  x > y we  have  (since  1 is  assumed 

1 1 


N 

odd)  that  0 (x)  = - JI  (e  - x)  = n (x  - e .)  JI  (x  - yj 

J-l  3 J-l  J £=1  3 


as  N -*■  °°. 


Case 


_2_>  - a,  i^  = 2m,  m £ 1.  Once  again  we  have  6^  < b 


and 


(a^,0^)  n I'  = 0.  Set 


N ®~1  m 

“ n [(y2j+1  + U)2j+1^  ~ X1  H ~ w2j^  ” *!•  Again  it  follows 


J-0 


J-l 


that  we  can  select  an  ^ such  that  Implies  that  0^^)  > 0, 

1 5 i s ip  01(x)  >0  for  x i (which  implies  that  0^(x)  > 0 for 
all  x £ X)  and  30^  £ p(l).  Setting  = y^  + oo^,  for  i odd  and 

- y±  - U±  for  i even,  we  have  for  y^,  1 S i s ip  that 


: 


f 
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l/M*(y  ) = l/NC-i)1"1^  n (n  - y,). 

j-1  J 


Again,  using  (2)  and  (3)  we  have 


(since  i^  is  even) 


. m _i  ^ y 

that  1/N0"(y  ) - f(y  ) n (y  - y ) n (y*  - yj 

j-ij+l  J j-1  J 1 


as  N -*•  <*>.  Also,  for  x e X with  x > y , we  have  that 

*1 


I?(X)  = n (n.  - x)  = n (x  - n.)  -*•  n (x  - y 
j=i  J j-i  J j-i  j 


) as  N -*■  «.  Note  that 


0^  has  the  same  limit,  as  N -*■  °°,  for  each  of  these  two  cases. 

Next  we  consider  the  case  where  > a.  In  this  case,  the  contri- 
N N 

bution  of  0^  is  identical  with  that  of  0^_  for  r = 2,...,v. 

Thus,  we  consider  the  construction  of  0f  for  ^ar«^r)  where 

r - l,...,v.  Here  we  must  consider  an  additional  four  cases.  For 
convenience  set  i.  = ir_^  and  k * iy  - if  ^ so  that 

<ar,Br)  n I = {y±  7±  } = {y1+1,...,yt+k}. 

r-1  r 


Case  3.  «r  > a,  6r  < b,  r = 1, . . . ,v;  k = 2m  + 1,  m > 0.  Note  that 
in  this  case  y|  = y£+k  e I'  and  p(r)  s 2m  + 2.  Set 


m 


m 


0r(x)  - (yA+k  x>^^[(y£+2j+1  + o^j+i)  " xI^1[(yJi+2j  ” “4+2J*  ~ x1- 


(4) 


Once  again,  there  exists  an  N such  that  for  H > N we  have 

r r 


1 


12 


L 


1) 

0%)  > 0, 

X 

< a 

r 

r 

il) 

^r^i)  > °’ 

i 

- 1 + 

lii) 

0”(x)  > 0, 

X 

* V 

(Note  that  ill)  follows  from  the  even  number  of  linear  factors  in  (4), 


N 

all  of  which  are  negative).  Therefore,  0^(x)  > 0 for  all  x e X. 


Setting  P±  = yl+±  + m£+i  for  i odd  and  P±  = for  i 

« 

L 

even,  we  have  for  x £ a 

r 

i 


0^(x)  = n (PA  - x)(y'  - x)  ->■  n (y  - x)  (y^  - x)  as  N -*■  ®. 

1=1  l j*sir_i+1  3 r 


Similarly  for  x S Br  (since  k is  odd) 


0^(x)  = n (p  - x)(y'  - x)  ->■  n (x  - y.)(x  - y'  ) 

i=!  Ct+K  j=ir_i+!  j r 


as  N . 


Finally,  for  y^,  i^_^  + 1 £ i £ i , we  have,  using  (2)  and  (3) 

(assuming  y^  = y^  for  some  0,  1 £ 6 £ r) , 


i-i  ,+l  1~ir-l~1 


1/N0^(y.)  - l/(-l)  r 1 0)  n (p,  - y.)  n (p.  - y.) (y!  - y.) 

j«l  3 J-i-i  X+1  3 1 r 


i-1 


n (y1  - yJ  n (y1  - yJ 

,l-i  1 3 1=i+l  3 1 


i-i  -1 
r 


5* 


fCyJ  n (y  - p ) it  (p  - y.) (y'  - y.) 
1 j-1  1 3 j-i-i  +1  3 1 9 1 


Y 

n (yi  - yt) 
J-i.  3 .. 
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.1‘r”1 

f(y4)  n (y,  - yJ  n 
J-l  J j-ir+l 


0-1 


(y,  - 


yt)  n (y^  - y ) n (y;  - y*) 
j-l  J j-e+i  J 6 


as  N -*■  00 . 


Case  4.  ar  > a,  < b,  k = 2m,  mil.  Note  that  in  this  case 
(ar»8r)  gives  no  contribution  to  I'.  Set 

K™  - ^o[(yW2J+1  + »t+2J+1>  - - »,+2J)  - xl  • 

As  before,  there  exists  an  Nf  such  that  N i N^  implies  0N(x)  > 0 
for  all  x c X.  Using  as  defined  in  Case  3,  we  have,  for  x ^ ot 

N k *r 

0 (x)  = n (p  - x)  -*■  n (y  - x)  as  N ■+  «,  and  for  x i 0 

j”l  J J-i^+1  J r 

N k *r 

0 (x)  ■ n (p  - x)  -*■  n (x  - y ) as  N ■+  Finally,  for  y . , 

j=i  J j-ir_i+l  J ± 

^.1  + 1 * 1 s ir.  we  have,  using  (2)  and  (3) . 


„ Vt  i-i  +1  1"ir-r1 

i/^0”(»)  * y*(-i)  r_1  0)  n (p 

j-l  J 


y.)  n 

j-i-i  .+i 

J r-l 


(P4  " 


yi> 


i-1 


n (yi  “ V n (yi  ■ V n lyi  * yj 

. i-i  J i-i+i  J 1 i-i  3 1 


^‘r-l-1 


f(y4)  n (y  - p ) n (p  - y ) 
j-l  1 • 3 j-i-i^ti  J 1 
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f(yj  n (y  - y ) n (y. 
J-l  J j-i  +1  J 


yj  n |y!  - y I as  N <». 
j-l  J 1 


Case  5.  > a,  0r  - b,  lr  - i ^ - k - 2m  + 1,  m 2 0.  Note  that 

In  this  case  r - v,  = p and  (av>6v)  gives  no  contribution  to  I'. 


Set  0v(x)  = ^[(yl+2j+1  + “4+2j+i>  - x]jf1[(yJi+2j  ' <u£+2j)  " x]  where 

is  defined  by  (3).  As  before,  there  exists  an  Nr  such  that  N 2 
N 

implies  0^(x)  0 for  all  x e X.  Using  as  defined  in  Case  3,  we 


for  x s a , 0 (x)  = n (p  - x)  -*■  n (y,  - x) 

^ J 1 j _ j .1  1 


as  N -*■  =». 


i=i  ,+l 
v-1 


Also,  for  y^  i ^ + 1 < i 


* N / i*i„-i+1  * , v-1  Y 

1^0  (y  ) - yV(- 1)  0).  n (p.  - y ) -*■  f(y  ) n (y.  - y ) n |y' 

j=l  3 1 1 j-l  1 J j=l  j 

A * 

as  N -*■  <»  using  the  same  argument  as  before. 


Case  6.  c»r  > a,  Br  = b,  (r  = v,  iy  - p) , p - ir_^  - k = 2m,  m 2 1. 
Again,  the  interval  (a  ,6  ) gives  no  contribution  to  I’.  Set 


Vx)  " j;o[(yi+2j+l  + “£+2j+l)  - x]j;i[(y£+2j  ■ “W  - x]-  **  before’ 

4 K 

we  have  0^(x)  > 0 for  all  x c X.  0N(x)  n (y  •-  x)  as  N -*■  », 

j-i  V 
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for  x s av.  For  y±,  + 1 * i * w. 


L N —1  Y 

lM  Cy.)  f(y.)  n (y  - y ) n |y!  - y |.  Recalling  that  p (s)  is 

j-1  J j-1  J 

the  number  of  zeros  of  q*(x)  in  the  interval  (a  ,0  ),  s - 

s s 

we  observe  that  in  Cases  4,  5 and  6,  30^(x)  ^ p(r).  Therefore,  we  set 

V N N 

pN(x)  = II  0 (x) , where  0 (x)  is  constructed  with  respect  to  the 
s=l  S 


K 


<«r 


interval  (a  ,6  ),  s = l,...,v  as  described  above  (depending  upon 
s s 

(a  »B  )).  Since  30N(s)  2 p(s),  these  intervals  (by  construction)  are 

S 8 S 

pairwise  disjoint  and  | | q*  | | - 1,  we  have  that  3p  < 3q*.  Also, 

A N 

there  exists  an  N such  that  N £ N*  implies  that  0N(x)  > 0 for  all 

s 

x e X,  s = l,...,v  so  that  pN  e(k)  for  N £ N*.  Furthermore,  for 
x c X ~ I, 


iim  „ , 

N-~  NPN(X) 


(iimi)  !! 

\N-x»  " / 8-l 


llm"N 

N-X»  0g(x) 


/ 1 1 

i Y i ,W 

“ flimi 

1 n - x 

n y!  - x - 0 

since  the 

Mi-i  1 

i-1  1 / 

quantity  in  parentheses  on 

the  right  is 

not  equal  to  zero. 

Also,  since 

IU  r'  dist(®,X  ~ I)  * c > 0,  we  have  that  the  above  convergence  is  uniform  in 
X ~ I.  Finally,  for  y^  e I (say  y^  e (°r>8r))  we  have  that 


— - 


i * l,...,p  and  < A holds  for  all  x e X ~ I we  have  that 

| |y  “ jj~|  lx  < A (recall  that  | £/xyl  - A on  X ~ I) , which  is  our 


desired  contradiction.  Hence,  we  must  have  c*  > 0.  Thus,  q*(x)  > 0 


c * 

V c ^ 

for  all  x e X.  Furthermore, ► — r uniformly  in  X as  v •+•  «■>  so 

\ q 


A This,  in  turn  implies  that 


p*  * q*/c*  is  our  desired  best  approximation  from  K. 


Finally,  we  would  like  to  close  this  section  with  results  on  characterization 
and  uniqueness.  Using  the  standard  argument  [17]  for  alternation  of  best  rational 
approximants  on  an  interval  we  have: 


i 

I . 

! I 


THEOREM  2.  Let  f e C ( X)  where  X is  a compact 

subset  of  the  real  line . Let  n ^ 0 be  a given  integer  and  set 

R®(X)  * (l/p(x):  p e n_,  p (x)  > 0,  for  all  x e X> . Assume  1/f  i R®(X) . 
n n — — — n 

Then  a_  necessary  and  sufficient  condition  that  1/p*  iis  .a  best  approxi- 
mation to  1/f  on  X from  R^(X)  is  that  the  error  curve 

e*(l/f)  ■ 1/f  - 1/p*  alternate  at  least  n + 1 times . 
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It  should  be  remarked  that  this  theorem  is  also  valid  for  Rm(X) 

n 

= (r  = p/q:  p e 11^,  q e JI^,  q > 0 on  X and  ||q  ||  = 1 1.  iitWTi  jin,  i Ui  i > in  , 

I Fijn-Tlir-u*  iji  il  JJm.iiiMi  i|l  in  uisr-ftm  in. 

On  the  question  of  uniqueness,  we  can  simply  quote  the  results  we 
desire  from  work  of  Cheney  and  Loeb  [6,  13].  Indeed,  they  actually 
prove  that  whenever  existence  holds  for  /(X)  then  so  does  uniqueness 
and  strong  uniqueness. 


THEOREM  3.  [16,  p.  94-95].  Let  f e C(X)  satisfy  f(x)  > 0 for 

all  x e X.  THen  there  exists  a unique  best  approximation  1/p*  to  1/f 
Rn(X>-  Furthermore,  there  exists  a positive  constant  y = y(f) 
such  that  for  each  1/p  e R^(X)  (strong  uniqueness) 


3.  COMPUTATION 

| 

Id  this  section,  we  wish  to  describe  three  possible  algorithms  for 
computing  best  approximations  to  1/f,  f(x)  > 0 for  all  x e X,  from 
R^(X) . The  algorithms  are  Remes,  differential  correction  and  a hybrid 
algorithm  which  is  a combination  of  the  first  two.  In  a future  paper 
we  shall  report  on  numerical  experiments  Involving  these  algorithms. 

In  what  follows,  we  shall  assume  X is  a finite  set. 

Remes  algorithm  has  been  widely  studied  and  appears  often  in  the 
literature.  Two  explicit  papers  where  the  Remes  algorithm  is  proposed 
for  calculating  best  rational  approximations  are  in  Cody,  Fraser  and 
Hart  [7],  and  Ralston  [15].  The  Remes  algorithm  consists  of  two  main 
operations : 

1)  the  solution  of  a nonlinear  system 
and  11)  the  exchange  of  a certain  set  of  points. 

It  is  known  that  in  general  the  nonlinear  system  may  have  many  solutions 
(and  sometimes  none  of  which  belong  to  R^(X))  [16,  p.  104],  [9].  Thus 
Remes  could  fail  to  run  due  to  its  inability  to  either  solve  this  system 
or  by  returning  a solution  to  this  system  which  is  not  in  the  class 
R^j(X)  . (A  second  problem  with  Remes  algorithm  will  be  mentioned  later. ) Even  if 
the  algorithm  is  able  to  solve  this  system  at  each  step  with  a solution  in  R^(X) , 
convergence  can  only  be  guaranteed  if  the  algorithm  was  initialized  with  a sufficiently 
"good"  starting  approximation  (for  the  case  X is  an  interval) . In  a recent  study  by  Lee 
and  Roberts  [12],  it  is  observed  that  the  Remes  algorithm  is  very  fast  when  it 
converges,  but  may  fail  to  converge.  Also,  observe  that  in  the  theory 


19 


we  are  considering  here  every  function  being  approximated  is  normal  In  the 
usual  sense. 

The  differential  correction  algorithm  was  originally  introduced  by 
Cheney  and  Loeb  [5].  This  algorithm  was  shown  to  have  very  desirable 
convergence  properties  by  Barrodale,  Powell  and  Roberts  in  [1]  and  a 
Fortran  listing  of  it  can  be  found  in  [11].  In  the  setting  considered 
here  this  particular  algorithm  possesses  guaranteed  (quadratic) 
convergence  to  the  desired  best  approximation  in  R^(X)  . However,  in 
practice,  if  X is  large  it  is  sometimes  necessary  to  solve  this 
problem  on  a subset  of  X and  then  initialize  the  full  problem  with 
this  solution.  This  is  due  to  the  fact  that  this  method  involves  a 
linear  programming  subroutine  which  is  sometimes  numerically  difficult 
to  solve  without  a good  initialization.  Also,  due  to  the  inclusion  of 
this  linear  programming  subroutine,  this  algorithm  is  quite  slow  (some 
19  times  slower  than  Remes  (when  Rentes  converges ) in  tests  done  in  [12]). 

The  final  algorithm  that  we  wish  to  mention  for  this  problem  is  a 
hybrid  of  the  above  two.  Precisely,  we  propose  to  replace  the  step  of 
the  Remes  algorithm  where  a nonlinear  system  is  solved  to  get  a best 
approximation  on  a ref erence  set  (smaller  than  X)  with  the  differential 
correction  algorithm  applied  to  this  reference  set  to  give  the  desired 
best  approximation  on  this  set.  This  method  will  eliminate  the  problems 
of  the  Remes  algorithm  associated  with  the  solution  of  the  nonlinear 
system  in  that  a best  approximation  on  the  reference  set  (which  is 
positive  on  the  reference  set)  will  be  found.  However,  there  is  still 


no  guarantee  that  the  best  approximation  on  a given  reference  set  found 
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by  the  differential  correction  algorithm  will  actually  belong  to  R^(X) 

(i.e.,  it  may  fail  to  be  positive  (or  defined)  on  some  points  of  X not 
in  the  reference  set).  In  fact,  we  have  encountered  such  examples  in 
testing  our  hybrid  algorithm  and  these  examples  have  given  rise  to  a second 
problem  in  the  Remes  algorithm.  Namely,  the  exchange  procedure  cycled. 

Thus,  it  is  also  necessary  to  modify  the  exchange  procedure.  At 
present  we  are  testing  two  modified  algorithms.  In  order  to  describe  these 
modifications,  let  us  assume  that  we  are  at  the  k-th  step  of  the  iteration 
and  suppose  that  the  best  approximation  r^  on  the  reference  set  (n  + 2 
points  from  X)  has  been  found  by  the  differential  correction  algorithm. 

The  first  modified  exchange  algorithm  is  as  follows.  Perform  a 
multiple  exchange  in  the  usual  manner  only  among  those  points  of  X where 
r^(x)  > 0 holds.  If  r^  is  not  the  best  approximation  on  the  set  of  points 
where  it  is  positive  then  a new  reference  set  is  obtained  and  the  algorithm 
proceeds  to  the  differential  correction  phase  to  find  the  best  approximation 


on  this  new  reference  set.  If  r^  is  the  best  approximation  on  the  set  of  points 
where  it  is  positive  then  terminate  the  algorithm  if  r^(x)  > 0 for  all  x e X 
(r^  is  the  desired  best  approximation  on  X)  or  adjoin  to  the  set  X^,  y e X 
where  r^(y)  = min{r^(x):  x e X}.  Note  that  r^ty)  <_  0 must  hold  in  this  case. 


Set  xk+i  = {y}  u X^  and  apply  the  differential  correction  algorithm  to  this 
set  of  n + 3 points  finding  the  best  approximation  r^+^  °n  it.  Next,  reduce 
X£+1  to  a subset  of  n + 2 points,  X^+^»  where  X^+^  is  chosen  so  that  alternation 
holds  on  X^+^.  Now  repeat  the  exchange  procedure  on  X^+^  with  respect  to  r^+^- 

The  second  modified  exchange  algorithm  is  basically  a reordering  of  the  above 


one.  In  particular,  if  r^  the  best  approximation  on  X^  is  positive  on  all  of  X 
then  we  proceed  with  a multiple  exchange  in  the  usual  manner.  If  r^  is  not 
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positive  on  all  of  X then  we  adjoin  y e X to  X^  precisely  as  above,  getting 
*£♦1  and  proceed  as  in  the  above  algorithm. 

One  can  prove  that  in  both  of  these  modified  algorithms,  cycling  cannot 
occur  and  that  global  convergence  holds  for  X finite  (i.e.,  error  of  approximation 
on  successive  reference  sets  strictly  increases).  In  a future  paper  we  shall 
give  a detailed  description  of  these  two  algorithms  and  report  on  the  results 
of  our  numerical  testing  of  them.  Also,  we  are  studying  the  extension  of  these 
ideas  to  R™(X)  and  will  also  report  on  this  at  that  time. 


The  second  author  wishes  to  acknowledge  a valuable  personal 
conversation  with  Professor  Dietrich  Braess  with  regard  to 
this  problem. 
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